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In the present paper we study the long wavelength and slow time scale behavior of a coasting 
beam in a resonator adopting a broad-band impedance model. Based on the renormalization group 
approach we derive a set of coupled evolution equations for the beam envelope distribution function 
and the resonator voltage amplitude. The equation for the resonator voltage amplitude is further 
transformed into a generalized Ginzburg-Landau equation. 

I. INTRODUCTION. 

So far nonlinear wave phenomena have received scant attention in the study of collective effects in charged particle 
beams. Considerable experimental and simulation data however exists suggesting that these phenomena should be 
included into the entire physical picture of beam propagation in accelerators and storage rings. 

A vast literature in the field of plasma physics is dedicated to the study of nonlinear wave-particle processes due to 
space charge interparticle forces. In high energy particle accelerators, where space charge forces are negligibly small, 
of particular interest is the coherent state of the beam under the influence of wakefields, or in the frequency domain, 
machine impedance. This state is highly nonlinear and depends on the interaction of nonlinear waves, involving some 
weak dissipative mechanisms that balance beam fluctuations driven be the wakefields. 

In a previous work Q, 0] we studied nonlinear behavior of a coasting beam under the influence of a resonator 
impedance. Starting from the gas-dynamic equations for longitudinal motion and using a renormalization group (RG) 
approach |s|, Q we found a set of coupled nonlinear equations for the beam density and resonator voltage. However, 
as is well-known, the hydrodynamic approximation is valid when the beam is close to a local equilibrium, which in a 
number of practically important cases may well be far from reality. 

The present paper, providing a complete kinetic description of the processes involved, is aimed to overcome the 
above mentioned difficulties. In what follows we study the longitudinal dynamics of a coasting beam in a resonator 
adopting a broad-band impedance model. We are interested in describing slow motion of beam patterns (droplets) 
neglecting fast oscillations of beam density and voltage on the resonator at a frequency close to the resonant frequency. 
We employ the RG method to derive amplitude equations governing the dynamics of slow processes. In Section II we 
obtain the desired equations for the longitudinal envelope distribution function and the amplitude of the resonator 
voltage. In Section III we proceed to transform the equation for the voltage amplitude into a generalized Ginzburg- 
Landau equation by solving explicitly the Vlasov equation for the envelope distribution function. Finally in Section 
IV we draw some conclusions resulting from the work performed. 

II. THE AMPLITUDE EQUATIONS. 

The starting point for the subsequent analysis is the system of equations: 
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for the longitudinal distribution function / (9, v; T) of an unbundled beam and the variation per turn of the voltage 
V (9; T) on a resonator. All dependent and independent variables, as well as free parameters in equations (|l|) are 
dimensionless and have been rescaled according to the relations: 

T = u s t ; w = — — = 1 + ; u= — ; 7 = — , (2) 
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Here uj s is the angular revolution frequency of the synchronous particle, AE is the energy error, ojr is the resonant 
frequency, Q is the quality factor of the resonator, R is the resonator shunt impedance and p is the uniform beam 
density distribution at the thermodynamic limit. Furthermore 

is the proportionality constant between the frequency deviation and the energy deviation of a non synchronous particle 
with respect to the synchronous one, while r\ — olm — 77 2 (&M - momentum compaction factor) is the phase slip 
coefficient. The voltage variation per turn V, the beam current / and the longitudinal distribution function / entering 
equations ([!]) have been rescaled as well from their actual values V a , I a and f a as follows: 

V a = 2eu s p -fRV ; I a = euo s p I ; /„ = p Q f (5) 

Let us introduce the Radon transform @], Q of the distribution function / (9, v; T) 

f {9, v- T) = ( d^F (9, Z;T)6[v-U {9, £ T)\ . (6) 



In the definition (^|) £ can be viewed as a Lagrange variable which is usually determined from the condition that the 
distribution function / (9,v;T) be equal to a specified distribution, say the equilibrium distribution for instance: 

f(6,v;T) = f (0 v = U(6,£;T). 

Substitution of eq. (||), into the system (Q) yields: 
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where the fact that the azimuth 9 is a slow variable (the dependence of F, U and V on 9 is through a stretched 
variable £ = e9) has been taken into account. Note that the system (0) resembles the set of gas-dynamic equations, 
governing the longitudinal motion of the beam. It bears however, additional information about the velocity distribu- 
tion, embedded in the dependence on the Lagrange variable £, and takes into account its overall effect through the 
integral on the right hand side of the third equation. 

We next examine the solution of the system of equations (0) order by order in the formal small parameter e by 
carrying out a naive perturbation expansion. The zero order solution (stationary solution) is readily found to be: 

F =F (O ; U = U (O ; V = 

and in particular one can choose 

M0 = /o(0 ; u (£) = i + £. 

Combining the first order equations 
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Solving the first order equations one easily obtains: 

Vi (0; T)=E (9; T ) e 4WlAT + c.c. 

C/x (0, C; T) = Uo (0, £; T ) + A ^ (g;To) e^ AT + c.c. (8) 
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where 
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and the amplitudes T (0; T ) , u (0, £; T ) , i? (6*, £; T ) are yet unknown functions of 0, £ and the initial instant of 
time To. Proceeding further with the expansion in the formal parameter e we write down the second order equations 
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and by elimination of U 2 and F 2 from the third equation we obtain: 
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Solving the above equation and subsequently the two other equations for U 2 and F 2 we find the second order solution 
as follows: 
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In a way similar to the above we write the third order equations as 
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Solving the equation for V 3 
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that can be obtained by combining the third order equations, and subsequently solving the two other equations for 
U 3 and F 3 we obtain the third order solution: 
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Next we collect the secular terms that would contribute to the amplitude equations when applying the RG procedure. 
Setting now e = 1 we write down the part of the solution of the system (Q) that has to be renormalized 

F RG (9, £; T, T ) = F (0, £; T ) - AT^ [f (0, £ T ) U (0, £; T ) 
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Following Kunihiro H we represent the solution (12) as a family of trajectories or curves {3?t } = 
[Frg (To) , Urg (To) , Vrg (To)] , being parameterized with Tq. The RG equations are defined as the envelope equations 
for the one-parameter family {sftjb} 1 
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From the above definition ([l3|) it is straightforward to obtain the desired RG equations: 
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where 



d£ (T - T ) 



E(9;T) = E{6;T)e—< T . 

The final step consists in defining the envelope distribution function G (9, v; T) by the Radon transform 

G(9,v;T) = f dZF(9,t;T)s\v-U(9,Z;T) 



(15) 
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By virtue of (|15|) the system of RG equations ([L4|) is equivalent to the following system of equations for the envelope 
distribution function G(9,v;T) and the resonator voltage amplitude E (6;T): 
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The system of equations (|16|) and ( |17| ) provides a complete description of nonlinear particle-wave interaction. It 
governs slow processes of beam pattern dynamics through the evolution of the amplitude functions. In (16) one can 



immediately recognize the Vlasov equation for the envelope distribution function G (9, v; T) with the ponderomotive 
force, due to fast oscillations at frequency close to the resonant frequency. It may be worth noting that the system 
( |l6| ) and ( |Tt[ ) intrinsically contains the nonlinear Landau damping mechanism, a fact that will become apparent from 
the treatment in the next Section. 



III. DERIVATION OF THE GENERALIZED GINZBURG-LANDAU EQUATION. 

In order to solve equation ([16]) we perform a Fourier transform and obtain: 
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and the Fourier transform of a generic function g (9; T) is defined as 
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g* (k;Sl) = [ g (-k;-n)\*. 
Solving equation (Hq) perturbatively we represent its solution in the form: 



G(x) = (27r) 2 fo(v)5(x) + G( X ) ; G( X ) = ^G„(x) 
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where 
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The Fourier transform of equation ( [17j ) yields the linear dispersion relation 
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entering equation (17) can be computed in a strai ghtf orward manner. Substituting the solution (JT^)-(^TJ) into the 
above equations with the linear dispersion relation ( |22|) in hand, up to second order in G (k, v; f2), we find 



I (6;T) = 1-W E(9;T) 
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is the Landau damping factor. In the above calculations the equilibrium distribution function has been taken to be 
the Gaussian one 
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and <te is the r.m.s. of the energy error, proportional to the longitudinal beam temperature. By substitution of the 
expressions ( ^3|) and (|24| ) into equation (|l7j ) we arrive at the generalized Ginzburg-Landau equation: 
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It can be further cast to a simpler form by introducing the rescaled independent and dependent variables according 
to 

AT cj ( XT \ \X\E 
t=— ■ x= , . =19-T+ ; * 
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Then the generalized Ginzburg-Landau equation for the amplitude of the resonator voltage takes its final form: 

i- + ^ = - (l - -1 j _ + a7 _ - (1 - W ) 

- 061 (7 + *u ) - ax (7 + *^ ) (1 - 6xW) A (^) + mi6 2 lbaJ °% ^, (27) 



where 
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61 = ^ > 01 = a(l - 46i) , 
and the function VL (l^l) is given now by the simple expression: 

Wai*l) = (l + *7L)(l-^l*| 2 ) l*| 2 - (28) 

The generalized Ginzburg-Landau equation (p7| ) is known [Q to provide the basic framework for the study of 
many properties of non equilibrium systems, such as existence and interaction of coherent structures, generic onset of 
travelling wave disturbance in continuous media, appearance of chaos. Recent experimental and numerical evidence 
(see e.g. |2| and the references therein) shows that similar behavior is consistent with the propagation of charged 
particle beams, and the generalized Ginzburg-Landau equation we have derived could represent the appropriate 
analytical model to study the above mentioned phenomena. 



IV. CONCLUDING REMARKS. 



As a result of the investigation performed we have shown that a coasting beam under the influence of a resonator 
impedance exhibits spatial-temporal patterns modulated by an envelope (amplitude) function, which varies slowly 
compared to the fast time and short wavelength scales of the pattern itself. Extracting long wavelength and slow time 
scale behavior of the system we have derived a set of coupled nonlinear evolution equations for the beam envelope 
distribution function and voltage amplitude. We have further shown that the amplitude of the nonlinear wave satisfies 
a generalized Ginzburg-Landau equation. 

It is worthwhile to mention that the analytical framework presented here bears rather general features. It pro- 
vides complete kinetic description of slow, fully nonlinear particle-wave interaction process, and allows higher order 
corrections to the generalized Ginzburg-Landau equation to be taken into account. 
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